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The Surface Green Function Matching analysis (SGFM) is used to study the normal modes of the 
, , interface oscillations between two non-mixed fluids by considering the difference in their densities 

and viscosities. The limiting case of viscous-inviscid system is used for comparison. The role of the 
viscosity and the density ratios on the momentum exchange and on the polarization of the surface 
modes is analyzed. 
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I. INTRODUCTION 

i 

The theory of surface waves in fluids is usually treated using the Orr-Sommerfeld equation obtained from potential 
method pL This procedure is useful to find the characteristics of the wave, such as dispersion relation and damping 
but becomes rather complicated when some other features are needed such as polarization and density of modes. 
On the other hand, the inclusion of the viscosities of all of the media make difficult to understand the physics of the 
^»- > , interface. Broadly speaking, the dual effects of viscosity is well known [Q: to dissipate the energy of any disturbance, 
but also it has the more complicated effect of diffusing momentum. At present, the theory for viscous cases is not 
nearly as complete or general as for inviscid cases and it provides only a partial understanding of the role of viscosity 
in such systems. 

y—( ' A suitable formalism for including all the viscosities with great ease in non-homogeneous systems studying the 
. response function has been developed elsewhere [|]-||] . This formalism, the method of Surface Green Function Matching 
lO ' (SGFM), has been extensively used to study various inhomogeneous problems involving surface waves at solid surfaces, 
both free solid surface (interface between vacuum and solid) and solid-solid interfaces ||. It has also been used in 
interface involving fluids [^JQ] as is the case of solid-fluid interface and even fluid-fluid interface (this last case analyzed 
to give an unified treatment of waves in solids and fluids which seem to be apparently unconnected problems) . As 
far as we know, there are no previous works where the SGFM have been applied to the hydrodynamics problems as 
[ these authors suggested. 

The aim of this paper is to apply the SGFM to the study of the physical characteristics (dispersion relation, damping 
and polarization) on the interface normal modes of two fluids at rest giving insight of the mechanisms of momentum 
exchange through the interface for different ratios between the viscosity and density of the two media. 

In the next section a brief outline of the main points of the SGFM is given for the fluid-fluid system at rest, 
highlighting the considerations made in the solution of the problem. Section III is devoted to the physical analysis 
Oh 1 of the polarization of the modes and the momentum exchange across the interface. In section IV it is carried out a 
numerical evaluation considering the physical interpretation of the terms and the results for pair of fluids which are 
analyzed as illustration. Finally some conclusions are outlined. 



II. SGFM FOR TWO VISCOUS FLUIDS INCLUDING SURFACE EFFECTS 

The formal development of the SGFM method has been fully explained elsewhere |B|1 and in particular the 
treatment of matching with discontinuities S, suitable for the case of two non mixed fluids where the interface has 
special effects not seen in the liquid bulk. Mathematical and formal details can be found elsewhere M and need not 
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be repeated here. It is only necessary to add that in fluid-fluid interfaces it is better to work with the velocity of the 
fluid particle in agreement with the Navier-Stokes equation, instead of the fluid particle deformation, suitable when 
solids are present. 

Consider a system formed by a fluid Mi for z < and a fluid Mi for z > 0, both of them at rest. It has a planar 
interface at z = 0. Analysing first each bulk media individually to prepare its description in a suitable way for the 
eventual matching at the interface, the coordinate system will be choosen considering the planes z =const as those of 
interest. The notation will be for coordinates r = (p, z), k = (n, q) where p and K are 2D vectors. 

As explained in ||, the SGFM start with the knowledge of the Green function G.F of the excitation studied in 
each bulk material constituent. Then, it is needed to analyze the physical model for the excitation to perform later 
matching at the interface. 

Now, to know the G.F of each bulk media, the 3D differential equations of hydrodynamics are the starting point. 
The fluids are usually treated as incompressible and described with the Navier-Stokes equation. However, as explained 
in H, it proves convenient here to give the theory for compressible fluids, even if compressibility effects are ultimately 
neglected. Then, the equation of mass for isoentropic processes, and the momentum conservation equation that 
govern the fluid motion are linearized by neglecting all nonlinear terms in disturbance quantities. They may be 
written, respectively, as 

~p(r,t) + pV-V(r,t)=0 (1) 

9 ^ V7„/„ ^ , „Yt2t/Y_ ^ , (.J , V 



Pg- t V(r, t) = -Vp(r, t) + ijV V(r, t) + [if + VV • V(r, t) (2) 

where c, p, p, rj and rf are the velocity of sound, dynamical pressure, equillibrium density, shear and bulk viscosities 
respectively, all of them considered as constants in each medium. V(r,t) is the velocity of the fluid. We neglected 
the external forces and supposed that the perturbation is small enough to neglect the convective term for pressure in 
(§• 

All space and time dependent quantities will be Fourier transformed according to exp[i(n ■ p — uit)] where lu is 
a frequency. Then, for surface wave propagation, the amplitudes are functions of (k, uj) on one hand and of z on 
the other. This z dependence is due to the fact that there is no spatial invariance in this direction and the Fourier 
transform can not be accomplished. Green functions, including the ones for the bulk material constituents, are then 
conveniently expressed as G(k,lu; z, z') or, simply, as G(z,z'), with (k,lo) understood everywhere. 

Time Fourier transform will be implied now on. From eq. (|l|) it is obtained p(r,oj) = (pc 2 /iw)V • V(r,uj) : which 
putted in eq. (||) gives rise to 

ipwVi(r,w) + (f - r]) — V ■ V{r,u) + r/V 2 ^(r,w) = (3) 

OXi 

with i = x,y, z and 

r = -^+(V + V) (4) 

IU! \ 3 / 

as the system of equations which couples the velocity components. This system must be solved as a whole as it can 
not be decoupled in the general case. 

The actual G(z, z') of each bulk media considered separately as infinity can be obtained in different ways but using, 
for instance, the Fourier transform 3D, it yields for the G.F ||: 



G(Jfe,w) ' 



ipuj — rjk 2 



(T-7])kk 



ipu — Tk 2 



(5) 



where I is the unit matrix and kk is a diadic product of the wave vector. 
Its poles 



/ • \ 1/2 / ■ \ 1/2 

qi = -n ) q t = -« J (6) 

describe the transverse and longitudinal modes of the infinite medium. In (|^) the incompressible fluid can be considered 
taking (T — » oo) and the proper limit is achieved. 



2 



There is no physical reason for the preference of a particular direction in the xy-plane. This spatial symmetry of 
the system allows us to define, for instance, k = (0, k, q) without loosing generality but getting simplification of the 
calculations. 

Note that qj — * i\n\ if the compressibility is neglected, see eq. (|]), given rise to a vanishing longitudinal mode. So, 
the qj pole describes the longitudinal mode due to the compressibility of the media. 

Let Gg be the Green function (G.F) of the surface system just defined and Qs its surface projection. Let Qg 1 
be the reciprocal of Qs in the two-dimensional p or k space. This is the central object in the SGFM analysis. In 
particular, knowing QZ it is possible to find the surface mode dispersion relation (SMDR) and the density of modes 
of the surface system B . It is important to stress that the secular equation for the SMDR, namely 

det Qs 1 =0 (7) 

expresses the continuity of the velocity and the stress components transmitted across z = 0. This is where the physics 
of the surface effects comes into the picture. These effects introduce changes in the stress components transmitted 
across the interface and are ultimately measured by some surface tensor ms whose physical meaning is that ms, 
acting on the velocity field V, yields the extra forces per unit area transmitted across the interface. 
Let us call Q$q to Qg 1 in the absence of such surface effects, then one finds @ 

g- s l = g s l + ms (8) 

Thus the problem is to find ms for the surface effects one wishes to study. It will be included in this case only the 
surface tension 7 according to Laplace's Law. It can be deduced p] that 
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(9) 



There is a little difference between the former expression and the expression obtained in |3j according to the fact 
that here the velocity of the fluid particle is considered instead the fluid particle deformation. 

Then, eq. (||) expresses the continuity of the velocities and the stress components transmitted across the interface at 
z = 0. Knowing Q^ 1 one can find the dispersion relation of the surface modes (SMDR) through the secular equation 

d) 



III. PHYSICS AND POLARIZATION OF THE SURFACE MODES 

The construction of Q~so is explained in ||. The result, after adding (||), is 

(10) 

where gg 1 is a 2 x 2 matrix and is the null vector 1x2. gg 1 has components 



QS 1 



mqti + ? 72qt2 
llgs 1 ! 



(gS 11 = ~2Z h_ 2 _ i 11 

K 2 +q n q tl K 2 + qi 2 q t2 

(g s 22 = — t + —_. — (12) 

« +quqti K +qi2qt2 iu} 

(Ss 1 )!* - -(gs 1 ) 21 = ( 2 ^ + 2in Vl ) - ( 2 P f K + 2i. m ) (13) 
\k 2 + quq t i / \k 2 + qi 2 qt2 / 

We shall refer to the modes as sagittal or S polarized with (0,Vy,Vz), transverse tangent or TT(Vx, 0, 0), longitudinal 
or L(0, Vy, 0) and transverse normal or TN(0, 0, Vz) modes, according to the component of the velocity they have. 

Now, on using ( |l0| ) in (f7|) the factorisation of the (^5 X )n matrix element yields a TT mode which does not interact 
with the others, whose dispersion relation is 

*7iqti + 7 ?2qt 2 = (14) 

and has x-axis polarization. 
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It is easily seen according to (Q) that the TT mode has no solution but as stressed in Q , it does contribute to the 
density of modes and therefore plays a non trivial role in the physical properties of the interface. This mode exists 
but it is not a stationary one if there is other surface effects considered || . 

The rest of ( |Io|) yields the secular equation 



detgg 1 =0 



(15) 



It gives a sagittal mode with polarization S(0, Vy t Vz) and surface tension included. It will be analyzed in the following. 

The factor (gg^ii, see eq. (0, represents the surface movement component in y direction due to compressibility 
of the media while (gg ^22 is a z direction surface movement. The factor (gg x )i2 represents a coupling between y 
and z movements giving rise to an S polarization mode. It means that the surface has both horizontal and vertical 
movements. In other words, the surface particles move in a kind of circular orbits depending of its phase difference. 

On the other hand, there are no important velocities in our system, then compressibility can be neglected as 
described in (|] and we will discuss whether the S polarization remains or not. Putting qj = i\n\ in (|Tl|), ( p^ ) and 
(O) it is obtained 



(SS 1 )! 



PlU>l\K\ 



p 2 UJl\K\ 



K 2 + i|K|q tl K 2 +i|At|q t2 

/ -In _ Pi^qti P2^q t2 

\Ss 22 — 2 i -I I 1 2 i -I I 

K +«Mqti « +»|«|qt2 



7«r 
i oj 



(gs 1 



s m 



-(gs 1 ) 2 ! 



!| K lqti 



2iK?7i — 



p 2 UJK 



i| K lqt2 



2iK?7i 



First of all let us consider the special case where the viscosity of Mi is neglected. If we put r\\ = in (|lf 
is obtained gti ~~ * 00 an d hence 



g S 



K 2 + i\n\q t2 

P 2UJK 2iKT,2 



P2^K 



*l«|qt2 



PIU 

i \k\ 



i l K lqt2 
P2^qt 2 



2iKT]2 



•*l«Rt2 



i w 



(16) 
(17) 
(18) 

(19) 



Note that even though the viscosity of one of the constitution! media was neglected, the coupled factor remains 
due to the nonzero viscosity of the other fluid. So, in this limit this mode remains as sagittal S(0, Vy, Vz) exhibiting 
movements in y- and z-axis for the surface particles. The fluid was taken as incompressible but there is still a 
component of velocity on y-direction. As far as we know, nobody has ever pointed out this fact clear, except Lucassen 
in his works J7|-p^|, who considered incompressible fluids, but the movement in j/-axis was due to active materials on 
the surface, no as an effect of viscosity. This coupling of movements could be responsible for a less wavelength and a 
bigger dissipation as it will be seen later. It is in this direction where the viscosity plays an important role. 

There is more information in eqs. (16)-([l8|). If viscosities are neglected (771 — rj 2 = 0) in (gg 1 ) 11 and (gg 1 ) 12 these 
expressions vanish because q^! — * 00, q^ 2 — > 00, but doing the same in (gg )22 does not yield a vanishing result. This 
leads to: 



gs 





p\u? + p 2 oJ 2 - 7« 2 |k| 



(20) 



The non zero component of Eq. (^0|) is a generalization of the Kelvin dispersion relation when the density of the 
upper medium is included. 

Then, it is concluded that when there is no viscosities the S polarized mode becomes TN mode because the coupled 
factor disappears and only remains (gg 1 )22- Indeed, if at least one of the viscosities is considered the coupled factor 
appears giving rise to the S mode, (see eq. jl9|)). Hence, the viscosity is the fundamental force which couples different 
modes among them. 

Furthermore, the viscosity is the main cause of momentum exchange between the two media through the surface 
on the y-direction movement. Note that the longitudinal component movement disappears according to (gg x )n — > 
when 771 and 7/2 are neglected. On the other hand, in ( p"5| ) the transverse normal movement described by (g s ~ 1 )22 
exists because of the densities and viscosities of the media, (see eq. (|l7|)). When the viscosities are neglected as in 
(pp|), the normal component movement still remains because of the densities of the fluids. Hence, in the case of the 
z-axis movement both the viscosities and densities are important for the exchange of momentum. 
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These results are in agreement to the fact that when the interface particle moves according to the longitudinal 
mode it remains on the plane 2 = and the viscosities are the only way for the two media to interact, but when the 
interface particle moves according to the transverse normal mode it goes into each medium sometimes at z > and 
other at z < and then the inertial effects of the media become important according to their densities. 

Expressions (|T5|)-(|T§|) also recover the Kelvin equation for an ideal fluid with free surface, (see references in 
Neglecting the viscosities, and setting p\ = it is obtained 

It can be seen that our formalism not only recovers the expression for the Kelvin equation, but also recovers the 
z-polarization of that mode. 

After this analysis one can return back to the problem for both viscous fluids. From (|il5|)-(|l8|) it is obtained the 
secular equation for the SMDR 

^ 2 [(Pi + P2)Oiqt 2 + PaQti) ~ MOi ~ P2) 2 ] + 1k 2 \k\ [pi(\k\ - q t2 ) + - Qti)] + 

+ 4k 2 |k|(?7 2 - 77i) 2 (|k| - q t i)(N - go) + Wwfe - t?i)(pi|k| - p 2 \n\ - Piq t2 + P2<hi) = (22) 

with the following new definition of q^ = (k 2 — ipuj/r/) 1 / 2 . This expression, which we recall corresponds to two 
viscous non mixed incompressible fluids, can also be accomplished applying the potential method, although using 
that formalism it is rather difficult to obtain the polarization of the modes. 

This is the equation to be used to study the modes if one includes both viscosities and surface tension effects for 
incompressible fluids. Expression (H|) was reported in || to study the surface waves at the interface between a solid 
and a fluid. They neglected the surface tension. One of the aim of this paper is to compare this theory with the 
theory which just take into account only one of the viscosities. From expression (|l9|) it is not difficult to achive the 
SMDR for the viscous-inviscid fluid interface 

- Lu 2 p 2 {p\ + P2) + 7« 2 |k|p2 + 4k 2 |k|77 2 (|k| - q t2 ) - Mp 2 uK 2 r\ 2 = (23) 

which reduces to equation (2.5) of H when 7 = and will be evaluated in the next section along with (^) for the 
viscous case. 





p 2 u? - 7k 2 |k| 



(21) 



IV. RESULTS OF THE NUMERICAL EVALUATION 



In order to make a numerical study the following quantities of length and time for nondimensionalization will be 
taken: 

time by Tq = ' 



P2I 2 



2 

length by L = — (24) 

P27 

The dispersion relation ( p2[ ) becomes 

u? [(1 + Q)(q^ + Qql 2 ) - |k|(1 - Q) 2 ] + « 2 |«| [|«|(1 + Q) - - Qq t2 ] + 
+4 K 2 |«|(1 - N) 2 (\k\ - qtiXM - <ff 2 ) + Ain 2 uj{\ - N) - Q) + cfo - Qql 2 ] = (25) 



for viscous fluids and eq. (|23| ) gives rise to 

-lu 2 (1 + Q) + k 2 \k\ + 4k 2 |k|(|k| - cft 2 ) - Aiun 2 = (26) 
for viscous-inviscid case, where Q = pi/p 2 , N = 171 / 772 and 

qu=^K 2 -i^y 2 (27) 
ql 2 = ( K 2 - tco) 1/2 (28) 
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Then the characteristics of the system will be studied by its SMDR with real values of the frequency u>. Let us 
allow k to be complex, its real part is 2ir times the inverse of the wavelength and the imaginary part is the distance 
damping coefficient (3 related with the viscosities of the media. The dimensionless parameters are k — 2ir/L and 
lu = 2it/T . 

Fig. m shows the SMDR for Q = 0.8. There is one mode which decreases its wavelength A and increases its distance 
damping coefficient [3 with increasing frequency at a fixed value of the parameter N . It is also shown that when the 
viscosity ratio N is increased the wavelength lightly decreases at any frequency. The curves split bigger at higher 
frequencies. On the other hand [3 increases with increasing N . Also it was plotted the curves obtained with N = 
from eq. (^6|) which means zero viscosity of the medium Mi . It can be seen that the theory which includes all the 
viscosities predicts small A and bigger (3 for a fixed u> with respect to the N = case. 

Fig. shows the dependence of K and (3 with respect to the variation of the density ratio Q at a fixed value of the 
frequency and viscosity ratio. It can be seen that A decreases when Q increases at a fixed N. This was deduced by 
Taylor in his study of the ripple formation on an infinitely thick viscous circular jet but neglecting the air viscosity. 
References are given in 0. We now prove that this is also true when both viscosities are considered. Also A decreases 
at a fixed Q when the viscosity ratio N takes higher values. So, the effect of the viscosity of medium M\ reinforces 
the effect produced by the density and it can be stated that the smaller wavelength will be obtained when Q and N 
are both bigger. It is also plotted the curve with N = corresponding to the viscous-inviscid case. It can be seen that 
the wavelength is always smaller in the case N (viscous- viscous case) . The curves split bigger as Q increases and 
(3 grows rapidly at low values of Q for a fixed value of N and tends to saturation for higher values of Q. This small 
variation of (3 with the variation of Q even at a fixed value of the viscosity ratio N reinforces the idea of the density 
as another mechanism of momentum exchange between the two media through the interface. It not only produces 
smaller wavelengths, but also produces lightly bigger distance damping coefficients (3. 

On the other hand, the distance damping coefficient (3 also increases at higher values of the viscosity ratio N at 
a fixed Q. It was also plotted the curves at TV = 0. It is interesting to note that the theory of viscous-inviscid case 
predicts a small decrease of the distance damping coefficient with increasing density ratio Q. This is in accordance 
to the fact that setting N = means to neglect the momentum exchange through the interface by the viscosity and 
raising Q represents to increase the dynamic properties of the surface given rise to a bigger distance for the wave to 
travel before vanishing. 

Fig. H shows the variation of A and (3 with respect to the viscosity ratio for a fixed value of u> at three values of Q. 
It shows that as the viscosity ratio increases, the wavelength reduces rapidly first and tends to a limiting value for 
N > 1. The curves start in the value of A corresponding to the viscous-inviscid case. Also, for a fixed value of N the 
wavelength decreases as Q increases, in correspondence with Fig. [|. For the coefficient (3 it is seen that it raises for 
increasing N. 

To illustrate this theory for real fluid combinations there will be used three pairs of fluids: air/water, water/aniline 
and water/mercury. The parameters of these fluids at room temperature are: 



Element Density Viscosity Surface 

Tension 

(kg/m 3 ) (rnPa s) (mN/m) 



air 


1.21 


0.018 




water 


998 


0.890 


71.99 


mercury 


13500 


1.526 


485.48 


aniline 


1022 


3.847 


42.12 



Then, for the system air/water it is Q = 0.0012 and N = 0.0202. In this case the SMDR is plotted in Fig. |. It 
can be seen that there is no difference of the wavelength reported by viscous-viscous and viscous-inviscid cases due 
to the small values of the density and viscosity ratios but there is a small increase of f3 for all frequencies when the 
air viscosity is considered. 

However, the operating conditions in many gas turbine combustors and liquid-propellant rocket engines are such 
that the density and viscosity ratios are higher. It could be so also in water-oil emulsions and other problems where 
the interface between two fluids plays an important role. Then, the SMDR for water/aniline (Q = 0.977, N = 0.231) 
and water/mercury (Q = 0.074, N — 0.583) systems were also plotted. In the first case (water/aniline) the densities 
are very similar but the viscosity of the aniline is much bigger than the water viscosity. In the case of water/mercury 
the viscosities are near one half one another, but the density of the mercury is much bigger the density of the water. 
The SMDR has been plotted in Fig. || and [| Note than in both cases the wavelength decreases and the distance 
damping coefficient increases in a visible way when the viscosity of the medium Mi (water in both cases) is considered. 
The difference is bigger at high values of the frequency where the viscosity effects become important. Then, one can 
conclude that the inclusion of the viscosities of all media produces a substantial decrease of the wavelength if the 
viscosity ratio are big enough. It gives rise also to a bigger distance damping coefficient for the wave. 
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V. CONCLUSIONS 



In the present paper the close relationship between the properties of low amplitude surface waves propagation with 
the viscosity and density ratios, in a system of two non-mixed incompressible fluids at rest has been set out. The 
SGFM method was used to accomplish the dispersion relation and the full study of wave propagation by varying 
different parameters of the media. 

It was shown that the viscosity is a fundamental parameter for the coupling of different modes. It gives rise to an 
S polarization mode with y and z components of the movement of the particles on the surface. Also it was seen that 
the viscosity is the main force in producing momentum exchange in the longitudinal mode, but for the transverse 
normal mode both the viscosity and the density ratios are important to the momentum exchange. 

When considering surface modes, it was shown that only one of them is allowed and its wavelength is smaller when 
considering the viscosity of both media for fixed values of the density ratio. Also it was seen a characteristic variation 
of the distance damping coefficient when the viscosity of all media are included. On the other hand the increasing 
of the density ratio also reduces the wavelength and produces a lightly increase of the distance damping coefficient, 
then this factor is also important in reducing the wavelength of the surface waves. 

In order to see more real situations, three pair of fluids were analyzed and the importance of taking in consideration 
all the viscosities was shown. 
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FIG. 1. Dispersion relation of the surface mode for Q = 0.8. The upper part gives the wavelength k/k and the lower part 
the distance damping coefficient (3/k 

FIG. 2. Relation between wavelength and distance coefficient with respect to the density ratio at a fixed frequency for 
different values of viscosity ratios. The case N = is the viscous-inviscid case 

FIG. 3. Relation between wavelength and distance coefficient with respect to the viscosity ratio at a fixed frequency for 
different values of density ratios 

FIG. 4. Dispersion relation of the surface mode for the air/water system. In the legend it is especified the Mi as left and 
M2 as right in the combination M1/M2, i. e., air /water in this figure. 



FIG. 5. Dispersion relation of the surface mode for the water/aniline system 
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FIG. 6. Dispersion relation of the surface mode for the water/mercury system 
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